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Abstract. We present the results of a detailed anal- 
ysis of the newly discovered 15' large radio galaxy 
WNB 0313+683. It has been discovered in the WENSS 
and NVSS radio surveys, and has been identified with an 
optical galaxy at a redshift of 0.0901 ± 0.0002. The hnear 
size of the radio structure is 2.0 Mpc, which places it in 
the class of Giant Radio Galaxies (GRGs). 
Radio observations have been carried out with the WSRT 
at 1.4 GHz and with the 92-cm broadband system, with 
the 100-m Effelsbcrg telescope at 10.45 GHz and with the 
VLA at 1.4 GHz and 5 GHz. At 10.45 GHz the core is 
the dominant source structure, emitting ^ 25% of the to- 
tal flux. It has an inverted spectrum with a spectral in- 
dex of a = +0.42 ± 0.03 (S^ cx j/") between 1.4 GHz and 
10.45 GHz, suggesting a very compact structure, although 
we can not rule out variability of the core luminosity. 
The Rotation Measure distribution has been mapped us- 
ing a new method and has been found to be very uniformly 
distributed over the source. It is therefore probably galac- 
tic in origin, with only a small contribution from the (sur- 
roundings of the) source itself. 

From the spectral index distribution we have determined 
an upper limit on the spectral age of 1.4 ± 0.1 x 10^ yrs. 
The particle density of the ambient medium, using ram- 
pressure equilibrium at the hotspots, is 1.6 x 10^^ cm^^^ 
for the southern lobe and <; 5.8 x 10^ cm^'^ for the north- 
ern lobe. An independent measure of the external den- 
sity has been determined using the amount of depolar- 
ization towards the southern lobe. This gives a density, 
averaged along the line of sight, which is a factor of 10 
higher (2.0 x lO^'^ cm^"^) than the density near the head 



of the lobes found using the ram-pressure arguments. This 
discrepancy might at least partly be the result of a con- 
tribution from internal depolarization, which we can not 
exclude on basis of our radio data. 

From spectroscopical observations of the host galaxy we 
find that the Ha emission line has a broad component, 
and that the extinction must be large with colour index 
E{B - V) = 0.98 ± 0.10 mag. Since the galactic lati- 
tude is +9? 8, the extinction is probably mostly galactic 
in origin. We further find that WNB 0313+683 has a very 
large optical emission line fiux with respect to its esti- 
mated jet power, when compared with the correlation be- 
tween these two properties found by Rawlings & Saunders 
(1991). We argue that this, together with the relatively 
high radio power and the inverted radio spectrum of the 
radio core, is suggestive of a new phase of radio activity 
in WNB 0313+683. 



Key words: Galaxies: active; Galaxies: Individual: 
WNB 0313+683; Radio continuum: galaxies; Intergalactic 
medium 



1. Introduction 

Giant Radio Galaxies (GRGs, e.g. Saripalli et al. 1986, 
Subrahmanyan et al. 1996) are the largest members of 
the radio galaxy population, with a (projected) linear size 
^ 1.0 Mpc^. At low redshift (z ^ 0.3) some thirty of these 
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large sources are known. The most extreme case is the 

radio source 3C236 at a redshift of ^ 0.09, which has a 
projected hnear size of 5.7 Mpc (WiUis et al. 1974, Strom 
& Willis 1980, Barthel et al. 1984). 

GRGs are interesting objects to study. First, at low red- 
shift their angular size is several arcminutes or larger, 
which allows detailed studies of the different components 
of their radio structure, such as their jets (e.g. NGC 6251, 
Perley et al. 1984) and their lobe emission (e.g. Mack et 
al. 1997a) using a variety of radio instruments. Secondly, 
because of the large physical size of GRGs, they have ex- 
panded well out of the denser central regions of the clusters 
that they reside in, into a much less dense Intergalactic 
Medium (IGM). This makes GRGs a powerful and unique 
tool to probe the low-density medium at a large distance 
from their host galaxies, which is inaccessible to current 
X-ray instruments. Thirdly, their active galactic nucleus 
(AGN) must be approaching the endstage of its radio ac- 
tive phase. GRGs can therefore provide important infor- 
mation on the properties of old AGN. Lastly, several of 
these large sources have been identified with quasars (e.g. 
4C 34.47, Jagers et al. 1982) or Broad Line Radio Galaxies 
(e.g. 0319-1-412, de Bruyn 1989). This makes them impor- 
tant test cases for the orientation dependent unification 
schemes for radio loud AGN (e.g. Barthel 1989). These 
state that quasars and radio galaxies are mostly alike, but 
that quasars have their radio jets oriented closer to the 
line of sight than radio galaxies. A natural consequence 
of this is that quasars are not expected to have large pro- 
jected linear sizes. 

There are now many indications that the ambient medium 
of the radio lobes of GRGs has a very low density (^ 
10~^ — 10^'' cm^''). Probably the strongest evidence is 
provided by observations of the Rotation Measure (RM) 
and the depolarization towards these sources. The mea- 
sured RMs are small, with values which are usually < 
20 rad m-2 (e.g. Klein et al. 1994, Strom & WiUis 1980). 
Depolarization occurs only at wavelengths above 10-20 cm 
(Willis & Strom 1978, Strom & Willis 1980, Jiigcrs 1986, 
Klein et al. 1996). This combination can only be explained 
satisfactorily by a very low density of thermal electrons 
along the line of sight. 

Further indications for a low density environment result 
from spectral age analyses of GRGs. Typical spectral ages 
for GRGs have been found to be 10^ — 10^ yrs (e.g. Mack 
ct al. 1998), which translates into expansion velocities of 
0.01c — 0.1c. Using the assumption of ram pressure equi- 
librium at the head of the jet (e.g. Miley 1980), external 
densities of ^ 10^® — 10~^ cm~^ are commonly found for 
GRGs. This method suffers from many assumptions that 
have to be made, such as equipartition between the energy 
density of the relativistic particles and the magnetic field, 
the filling factor of the radio lobes, the fraction of energy 
in heavy particles, the area of the bow shock, and so on. 
It is therefore often criticized (e.g. Eilek et al. 1997) and 
its result should be interpreted with care. 



Many clusters have large and bright X-ray haloes, often 
extending to distances of more than 1 Mpc from their cen- 
ters and containing cooling flows (e.g. Fabian 1994). X-ray 
studies of some GRGs (e.g. NGC 6251; Mack et al. 1997b), 
however, have found only weak thermal X-ray emission 
around the host galaxies. GRGs are therefore probably 
not inside rich clusters with dense cores, and the lobes are 
thus not likely to be found in dense environments. Subrah- 
manyan et al. (1996) studied a small sample of GRGs on 
the southern hemisphere. By studying the surface density 
of optical galaxies in the neighbourhood of the GRG host 
galaxies using the UKST plates, they conclude that they 
do not reside in rich clusters. 

GRGs have relatively low radio powers (Saripalli et al. 
1986; Subrahmanyan et al. 1996), usually around or be- 
low the luminosity which divides FRI and FRII type ra- 
dio galaxies (- 2 x lO^^ W Hz"! at 178 MHz; Fanaroff 
& Riley 1974). Because of their large sizes and relatively 
low radio power, the surface brightness of GRGs is low. 
This is why they are so difficult to find in most radio 
surveys. In recent years, the Westerbork Northern Sky 
Survey (WENSS, Rengelink et al. 1997) has mapped the 
sky above -1-28° declination at a frequency of 327 MHz. 
WENSS has a sensitivity of 18 mJy (5(t) and a beanisize 
of 54" X 54" cosec(5, with 8 the declination. At low fre- 
quency, the bulk of emission of radio galaxies originates in 
the extended radio lobes. This makes the WENSS ideally 
suited to find large, low surface brightness radio galax- 
ies such as GRGs. The first discovery of such an object 
(WNB 1626+5153) has been reported by Rottgering et 
al. (1996). This initiated a project aimed at finding and 
studying a large, uniformly selected, sample of giant radio 
sources from the WENSS. 

Here we report on the discovery and subsequent analy- 
sis of the radio source WNB 0313+683, which we have 
identified as an FRII-type radio galaxy with a linear size 
of 2.0 Mpc. It has several remarkable properties, among 
which there is a large flux asymmetry of the radio lobes 
and a prominent, inverted spectrum, radio core. In Sect. 
2, we will present the radio and optical data that we have 
collected on this object. Section 3 describes a first analysis 
of these data, including a new way to measure Rotation 
Measures. It also presents some of the derived physical 
properties. In Sect. 4 we derive the advance velocities of 
the hotspots and the age of the source from a spectral 
index analysis. Section 5 then discusses the observed de- 
polarization towards WNB 0313+683. A discussion on the 
properties of the radio core is given in Sect. 6. We argue 
that WNB 0313+683 may currently be in a new phase of 
radio activity. Finally, a summary and our conclusions are 
presented in Sect. 7. 

2. Observations 

In this section we describe the observational data we have 
obtained for the source WNB 0313+683. An overview of 
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Table 1. 



Non-survey radio observations of WNB 0313+683 



Telescope Date Freq. Bandwidth Int. time r.m.s. 



WSRT Aug. 6 1995 1435 MHz 5 X 10 MHz 32 mins. (8 X 4) 0.3 mjy 

WSRT Oct. 27, Nov. 11 1995 343 MHz 8x5 MHz 9 hrs. (6 + 3) 0.5, 1.4, 0.7 mJy (Stokes V, Q, U) 

Effelsberg 100-m Dec. 10-27 1995 10.45 GHz 300 MHz 1 min beam-^ 1.3 mJy (Stokes /), 0.37 mJy (Stokes Q, f) 

VLA B-conf. Nov. 19 1995 1425 MHz 2 x 25 MHz 7 mins. Combined with C-conf. data 

VLA C-conf. Feb. 19 1996 1425 MHz 2 x 25 MHz 10 mins. 0.13 mJy (NA-weighting) 

VLA C-conf. Feb. 19 1996 4860 MHz 2 x 50 MHz 10 mins. 0.04 mJy (NA-weighting) 



radio observations which are not related to pubhcly avail- 
able surveys (i.e. the WENSS and the NVSS) can be found 
in Tab. 1. 



2.1. Survey data: WENSS and NVSS 

The radio source WNB 0313+683 was noticed as a 

15' large FRII-type radio galaxy in the WENSS and 
in the NRAO VLA Sky Survey (NVSS; Condon et al. 
1993). Figure 9a shows the radio contours from a map 
of WNB 0313+683 from the WENSS survey. Clearly vis- 
ible are the two bright regions at the outer edges of the 
radio lobes, the bridge connecting these regions, and a 
bright central lateral component extending to the north- 
west. The brightness and narrowness of the southern lobe 
suggests the presence of a radio jet. The northern lobe 
shows a protrusion at the top. From the WENSS map 
alone it is not clear if this protrusion is part of the ra- 
dio galaxy (i.e. a hotspot), or due to an unrelated radio 
source. The total integrated flux density at 327 MHz is 
3.12 ± 0.15 Jy. Flux densities for the different source com- 
ponents are given in Tab. 2. 

The angular resolution of the NVSS survey is 45", com- 
parable to that of the WENSS survey. The typical noise 
in the produced maps is 0.45 mJy beam~^. Further, the 
NVSS provides maps of the Stokes parameters Q and U . 
Figure 1 shows the radio map from the NVSS survey, with 
vectors of the observed £J-field superposed. The southern 
lobe and the central 'bulge' are still well detected. The 
bridge to the northern lobe is very faint, which is a first 
indication of its steep spectrum. The polarization angles 
in the southern lobe are more or less constant in orienta- 
tion; only near the core and the bulge there is a gradual 
change in polarization angle. 

The weak sources at RA 03''14'°34^5, Dec. 68°18'10" and 
RA 03'^13'"13?6, Dec. 68°21'37" (B1950.0) are most prob- 
ably unrelated background sources. 



2.2. WSRT Observations 

2.2.1. WSRT 1.4 GHz snapshot observations 

To find an accurate position of the radio core to fa- 
cilitate identifying the host galaxy, we have observed 
WNB 0313+683 with the Westerbork Synthesis Radio 
Telescope (WSRT) at 1.4 GHz. We obtained a short ob- 
servation on August 6, 1995, consisting of 8 x 4 minutes of 
integration, separated by typically 1.5 hours. Only a short 
(2 min) observation was made of the primary calibrator 
3C286, which did not allow accurate (^ 10%) flux cali- 
bration. Furthermore, the sparse UV-coverage (basically 
8 spikes at regular hour-angle intervals) is not adequate 
to obtain a good radio map of all the sources structures. 
Therefore, we do not quote any flux values from these ob- 
servations. 

Calibration and mapping of the data were done using the 
NFRA NEWSTAR package. The resulting map is presented 
in Fig. 2. A radio core is clearly detected and it is un- 
resolved. Its B1950.0 position is 03'^13™37!90 ± 0?05 in 
right ascension and +68°18'34'.'08±0'.'2 in dccHnation. The 
southern hotspot region is resolved and elongated along 
the radio axis. A radio jet has not been detected in these 
observations. 



2.2.2. WSRT 92-cm broadband observations 

We observed WNB 0313+683 using the 92-cm broadband 

system of the WSRT. to map the extended regions with 
low surface brightness in more detail. The 92-cm broad- 
band system uses 8 channels of 5 MHz, positioned between 
319 MHz (94.0 cm wavelength) and 380 MHz (78.9 cm 
wavelength). We observed WNB 0313+683 for 9 hours, 
distributed over two periods of 12 hours. The primary cal- 
ibrators 3C 48 and 3C 286 were observed, which we used 
for gain and phase calibration. We use the scale of Baars et 
al. (1977). Due to man-made radio interference, only 5 of 
the 8 channels (the ones between 325 AfHz and 360 MHz) 
could be used for mapping and further analysis of the 
data. Only one step of phase selfcalibration was necessary 
to improve the data quality to nearly the theoretical limit. 
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Fig. 1. Contour plots of the radio source WNB 0313+683 at 1400 MHz, from the NVSS survey, a Contours of total intensity 
and the i5-vcctors superposed. Contours are at -1.5,1.5,3,6,12,24,48,96,192 mjy beani"^. The length of the vectors represents 
the polarized intensity (1" = 0.2 mJy beam~^). b Contours of polarized intensity with the fractional polarization in greyscale. 
The greyscale ranges from 0% (white) to 50% (black). Contourlevels are at 1.2,2.4,4.8,9.6,19.2 mJy beam~^. 




Fig. 2. a Contour map of the 1.4 GHz WSRT snapshot observations. Contour levels are at -1,1,2,4,8,16,32,64 mJy beam" . 
b Overlay of the 5-GHz VLA radio map of the core of WNB 0313+683 (contours) and an optical image extracted from the 
Digitized Sky Survey (greyscale). The radio core is identified with a faint and extended galaxy-like object. 
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Due to possible confusion by weak background sources in 

the total intensity map, the best estimate of the noise is 
obtained from the Stokes V map. We measured a r.m.s. 
noise of ~ 1 mJy in each channel, and ^ 0.5 mJy in the 
average of the five useful channels. 

For each channel, maps were made in Stokes /, Q, and 
U parameters. The total power maps have been averaged 
into one single map, which is represented in Fig. 3. The 
sensitivity of the new observations is 6 times better than 
that of the WENSS survey. We detect additional emission 
outside of the lowest WENSS contours. 
Even at a wavelength of 92 cm there is still a substantial 
amount of polarization measured. We have mapped the 
polarized intensity in each individual channel, and aver- 
aged the results from the five good channels in one sin- 
gle map. The result is also shown in Fig. 3. Around the 
southern hotspot there are some residual sidelobes visible 
in the polarized intensity map. We decided not to remove 
them, because they fall mainly outside the total intensity 
contours and are therefore easily discernible. We have not 
plotted the polarization vectors, because the amount of ro- 
tation that has occured at such a low frequency deprives 
this of any physical meaning. 

2.3. Effelsberg observations 

To study the high-frequency morphology and polariza- 
tion properties of WNB 0313-1-683, it was observed at 
10.45 GHz with the 100-m Effelsberg telescope in Decem- 
ber 1996. At this frequency the beamsize is 69", which 
is close to the WENSS beam. Furthermore, since the ob- 
servations are done with a single-dish telescope, they are 
sensitive to large-scale structures, in contrary to most 
interferometers at this frequency. The observational and 
data reduction procedures are detailed by Gregorini et al. 
(1992). The data have been calibrated according to the 
scale of Baars et al. (1977). In addition, the maps have 
been CLEANcd, applying the algorithm described by Klein 
& Mack (1995). The final maps are shown in Fig. 4. The 
noise is 1.3 mJy beam~^ for total, and 0.37 mJy beam— 1 
for polarized intensity. The total intensity map is domi- 
nated by the core and the southern lobe. We measure a 
flux density of the core of 55 ± 3 mJy at 10.45 GHz. 
Polarized emission is detected in the radio lobes, but not 
towards the core. The southern radio lobe shows polariza- 
tion degrees of up to 40%. At this high frequency, the an- 
gles of the observed i?- vectors are very close to their true 
direction^, so that we do not have to correct for Faraday 
rotation to derive the direction of the projected magnetic 
field in this source. The observed £'-vc;(;tors are perpen- 
dicular to the radio axis, which means that the projected 
magnetic fleld must be largely parallel to the radio axis. 

^ Faraday rotation rotates the polarization angle over Ai9 = 
RM-\^ rad. For \RM\ < 50 rad m"^, A?? < 2°, which is below 
the accuracy of our observations. 



This is comparable to what is generally found in powerful 
FRII radio sources (e.g. Saikia & Salter 1988). 

2.4- VLA observations 

The source WNB 0313-^683 was observed with the VLA 
in its B and C configuration on November 19th 1995 
and February 19th 1996, respectively. We observed 
WNB 0313-1-683 for 7 minutes in the B configuration, and 
for 10 minutes in the C configuration, using two stan- 
dard 25-MHz IF bandwidths at frequencies of 1385 and 
1465 MHz. Furthermore, we made a 10-min observation in 
the C configuration with two 50- MHz IFs at frequencies of 
4835 and 4885 MHz. The phases were calibrated observ- 
ing the nearby source 0217-1-738, and the radio sources 
3C 286 and 3C 48 were used as primary flux density cali- 
brators. Calibration and mapping of the data was carried 
out with the NRAO AiPS package. To produce the final 
maps, the 1.4-GHz data from the B and C configurations 
were combined. The data have been mapped and selfcal- 
ibrated to solve for phase and gain variations during the 
observations. Total intensity and polarization maps have 
been produced, after CLEANing of the source structure us- 
ing Natural Weighting of the UV-data (see Figs. 5, south- 
ern lobe, and 6, northern lobe). At 1.4 GHz, both lobes 
are mapped in high detail. The southern lobe appears well 
confined; the width changes only by some 20% along its 
length, and it has sharp outer edges. The northern lobe 
is more diffuse and its hotspot is much less bright. The 
5-GHz observations only reveal an unresolved core and no 
other source structures. The flux density of the core at 
this frequency is 40 ± 3 mJy. 

Also shown in Figs. 5 and 6 are maps of the fractional 
polarization at 1.4 GHz. In the southern lobe, near the ra- 
dio core, there is a region of strong fractional polarization 
with values reaching 60-70%, which is close to the theo- 
retical maximum. This must therefore be a region with a 
highly ordered magnetic field. The orientation of the ob- 
served -E-field is here perpendicular to that in the south- 
ern lobe. Halfway the southern lobe there are two regions 
which also have a high fractional polarization. Although 
the outer edges of a radio lobe are often spuriously highly 
polarized, the extent of these two areas is too large to be 
ascribed to this effect, and we believe the increase in frac- 
tional polarization to be real. In the northern lobe only 
very little polarized structure is seen. The radio core is 
unpolarized at our detection limit. Its fractional polariza- 
tion is therefore below 5%, which is expected in the cores 
of radio galaxies (Saikia & Salter 1988). 



2.5. Optical Observations 

The host galaxy was first identified by making an over- 
lay of the WSRT 1.4-GHz map and the optical field ob- 
tained from the Digitized Sky Survey (the red POSS I 



6 



A.P. Schoenmakers et al.: WNB 0313+683: Analysis of a Newly Discovered Giant Radio Galaxy 




L 

00 13 30 00 

RIGHT ASCENSION (B1950) 



68 26 



^ 22 



O 
5 



o 
111 
Q 



25% 



20% 



15% 




10% 



03 15 00 14 30 



00 13 30 00 

RIGHT ASCENSION (B1950) 



12 30 



Fig. 3. a Contour map of the 92-cm broadband observations. The map shown here is the average of the 5 good channels in these 
observations, and represents the source at 343 MHz. Contour levels are at -4,4,8,16,32,64,128,256,512 mjy beam~^. b Contour 

map of the polarized intensity, averaged over the 5 good channels of the 92-cm broadband observations. Contour levels arc at 
5,7.5,10,15,20,30,50,100 mJy beam~^. The greyscale represents the fractional polarization, on a scale from 0% (white) to 25% 
(black). 



plates). At the position of the radio core a weak optical 
galaxy is seen (see Fig. 2). Its B1950.0 coordinates are 
03^13^38^33 ± 0?08 in RA and 68°18'34'.'3 ± I'.'O in dec- 
lination. The higher resohition 5-GHz VLA observations 
later confirmed the identification (see Fig. 2). 
We observed this galaxy with the 2.5-m INT telescope on 
La Palma, using the IDS spectrograph equipped with a 
Ikxlk TEK chip. A pixel scale of 3.12 A/pixel in wave- 
length was used, resulting in a total wavelength coverage 
of ^ 3100A. The spatial scale along the slit is O'.'84/pixel. 
The orientation of the slit was parallel to the radio axis, 
and a slitwidth of 2" was used. This gives a spectral reso- 
lution of loA. 

One 1200-s exposure was made in August 1995, with a 
central wavelength of 5500A. In October 1996, two addi- 
tional 600-s exposures were made at a central wavelength 
of 6OOOA. After the galaxy exposure, an exposure of a 
nearby F8-type star was made. This we used to correct 
for the atmospheric absorption bands. All spectra were 
reduced in the standard way using the longslit pack- 
age in the NOAO iraf data reduction software. Wave- 
length calibration was done using arc-lamp exposures and 
checked against the skylines on the original frames. Flux 
calibration was achieved by observations of several spec- 
trophotomctrical standard stars during the night. In the 
final extraction a 4" spatial aperture was used. This corre- 



sponds to 9 kpc at the rcdshift of the galaxy. The resulting 
spectrum is shown in Fig. 7. It shows a wealth of emis- 
sion lines, typical of a narrow-line active galactic nucleus 
(AGN), on top of a weak stellar continuum. Using the 
emsiion-lines we measure a redshift of 0.0901 ± 0.0002. 

3. Results 

In this section we determine several source characteris- 
tics from the available data. First, we discuss the over- 
all radio morphology and some properties of the radio 
core. Secondly, we produce spectral index maps using the 
WENSS, NVSS and Effelsberg data. Thirdly, we investi- 
gate the radio polarization properties of the source. We 
apply a new method for finding Rotation Measures and 
we present maps of the observed depolarization towards 
WNB 0313-1-683. Then, we derive some general physical 
properties of the source from the radio data, such as in- 
ternal energy densities and pressures. Lastly, we analyse 
the optical spectrum in more detail. 

3.1. Morphology, fluxes, and the radio core 

The morphology, although typical of an FRII radio galaxy, 
shows some interesting features especially at 327 MHz. As 
mentioned, the most prominent feature is the bright nar- 
row southern lobe, and the bright region at the end of 
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Fig. 4. a Contour map of the 10.45-GHz Effelsberg observations of WNB 0313+683, with the measured vectors superposed. 
Contours are at -4,4,8,16,32 mJy beam^^. The length of the vectors represents the polarized intensity (1" — 0.1 mjy beam~^). 
b Contour map of the polarized intensity at 10.45 GHz. Contours are at 1,2,4,8 mJy beam~^. The grayscale represents the 
fractional polarization and ranges from 0% (white) to 50% (black). 



it. The VLA observations have resolved this hotspot into 
two distinct bright compact components with a narrow 
spur of emission between them. It might be the case that 
WNB 0313+683 has two hotspots in its southern lobe, sim- 
ilar to what is seen, for instance, in the well studied case of 
Cygnus A (e.g. Carilli et al. 1991). The difference is that 
in WNB 0313+683 the two hotspots are well aligned with 
the radio axis. 

The northern lobe can be roughly divided into two parts: 
the brighter northernmost region and the much weaker 
extended emission that fills up the area between the core 
and the brighter part (the bridge). Although there is a 
hotspot, it is not as powerful as its southern counterpart. 
The radio lobe as a whole has a much more relaxed mor- 
phology. 

An important parameter to describe the structural evo- 
lution of radio galaxies is the ratio between the length 
and the width of the radio source, the so-called axial ratio 
(e.g. Leahy & Williams 1984). A low axial ratio indicates 
a rapid lateral expansion of the radio emitting plasma, 
most likely caused by a high internal pressure in the radio 
lobes, as compared to the pressure in the ambient medium 
(e.g. Begelmann & Cioffi 1989; Nath 1995). The axial ratio 
can best be measured at a low frequency, because of the 
increased visibility of the old, extended regions. Using the 
method described by Leahy & Williams (1984), we find a 



value of ^ 6.7 from our 92-cm WSRT observations. This is 
similar to the axial ratios of other GRGs (Subrahmanyan 
et al. 1996) which are typically between 4 and 11. It sug- 
gests that WNB 0313+683 has a similar internal pressure 
as other GRGs. It will be shown later that the equipar- 
tition pressures in the radio lobes of WNB 0313+683 are 
~ 10^^^ dyn cm^^, which is indeed within the range found 
for other GRGs (Subrahmanyan & Saripalli 1993, Mack et 
al. 1998). 

Table 2 gives the flux densities of the different components 
of WNB 0313+683 at different frequencies, from our own 
data and from the literature and on-line databases. Where 
possible, we have made a distinction between the north- 
ern lobe, the southern lobe and the core. If no value for 
the core is given, it is included in the flux density of the 
southern lobe. 

The NVSS flux density, although originating from in- 
terferometer data, agrees well with the value from the 
single-dish NRAO 1400-MHz survey (Condon & Brod- 
erick 1985). The beam of the latter is ~ 13' x 11' so 
that the three background sources visible in the NVSS 
map are incorporated in the source. We have estimated 
the flux density of WNB 0313+683 by subtracting the 
NVSS flux densities of the background sources from the 
NRAO 1400-MHz flux density. The NVSS flux density of 
WNB 0313+683 is still somewhat lower than this estimate, 
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Fig. 5. a Contour map of the southern radiolobe from our 1.4-GHz VLA observations. Contour levels are at 
-0.4,0.4,0.8,1.6,3.2,6.4,12.8,25.6,50.2 mjy beam"'^. Superposed are the electric vectors of the polarized emission. The length 

of the vectors is proportional to the polarized intensity (1" — 0.18 mJy beam~^). b Contour map of the polarized intensity. 
Contour levels are at 2.4,4.8,9.6,19.2,38.4 mJy beam~^. The greyscale represents the fractional polarization. The range is 0% 
(white) to 70% (black). 



which can be expected for a source of this angular extent. 
Figure 8 shows a plot of the values given in Tab. 2. 
Note the upturn of the spectrum of the northern lobe at 
10.4 GHz. A similar effect has been observed in the GRG 
1358+305 (Parma et al. 1996), which was also observed 
with the Effelsberg telescope at the same frequency. 
From Fig. 8, it appears that the radio core has an inverted 
spectrum. But since the data have been taken at different 
epochs, it can also be variable. To investigate this we use 
the 1.4-GHz WSRT and VLA observations, which were 
done some 3 — 6 months apart. Although the flux density 
calibration of the short 1.4-GHz WSRT observations was 
not very accurate, we compared the WSRT fluxes of three 
background sources in the field and of the core with the 
VLA observations. In case of no variability, the ratio of the 
VLA and WSRT flux densities should be the same within 
the errors. This is indeed the case for the three background 
sources, but the core flux ratio differs from the mean of 
the background sources by ^ 2a. To calculate the errors in 
these values we assumed a 2% uncertainty in the flux den- 
sity calibration of the VLA data and a 10% uncertainty in 
the WSRT data. Part of the deviation of the radio core flux 
ratio can be explained by the VLA observations picking 
up more of the extended emission surrounding the core. In 
the WSRT data, the UV-plane is poorly sampled and the 



map misses much of the extended structure. Therefore, 
we do not think that the found variability is significant, 
although we cannot rule it out either. Additional data is 
needed to put more stringent constraints on this. 
Neglecting the possible variability, we find a core spectral 
index of +0.42 ± 0.03 between 1400 and 10450 MHz from 
a least-square fit through the three available flux density 
measurements. The turnover in the spectrum of the core 
must be at a frequency above 10 GHz. 
We can estimate its size by assuming that the inverted 
spectrum is caused by synchrotron self-absorption. The 
equipartition angular size i^eg of the core is the size for 
which a source with given peak flux density and peak 
frequency is in equipartition (Scott & Readhead 1977). 
Assuming that the spectral index above the turnover is 
^ —0.75, we find -deq ^ lO^'"'. This result varies very 
little on the chosen spectral index. At the redshift of this 
source, 10""^" translates into a linear size of only 0.2 pc. 
The fraction of the total emission that comes from 
the core, or the core-ratio C (e.g. Orr & Browne 
1982), is a strongly increasing function of frequency for 
WNB 0313+683: At 1400 MHz, C = 0.03, at 4850 MHz, 
C = 0.12 and at 10450 MHz, C = 0.25. The median value 
for 3CR radio galaxies is only 0.002 at 8 GHz (Saikia & 
Kulkarni 1994). Values as high as --^ 0.2 at 8 GHz are nor- 




Fig. 6. a Contour plot of the northern radiolobe from our 1.4 GHz VLA observations. Contour levels are at -0.4,0.4,0.8,1.6,3.2,6.4 

mjy beam~^. Superposed are the electric vectors of the polarized emission. The length of the vectors is proportional to the 
polarized intensity (1" = 0.09 mJy beam^^). b Contour plot of the polarized intensity. Contour levels are at 2.4,3.6,4.8,7.2 
mJy beam~^. The greyscale represents the fractional polarization. The range is 0% (white) to 70% (black). 
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Fig. 7. The optical spectrum of the host 
galaxy of WNB 0313+683, corrected for at- 
mospheric absorption features. Details on 
the observation can be found in the text. 
The most prominent emission lines have 
been indicated. 



mally only seen in SCR quasars. Saripalli et al. (1997) find 

that in GRGs a typical value for C at 1.7 GHz is C ^ 0.01, 
and probably even somewhat larger at higher frequencies 
due to the usually steep spectrum of the radio lobes. Still, 
a core-ratio of ^ 0.2 at 8 GHz is very high, even for GRGs. 
The radio core prominence can be used as an orientation 
indicator in the light of relativistic beaming and unified 



scheme models (e.g. Giovannini et al. 1994, Morganti et 

al. 1997). Comparing the measured core radio power with 
the expected intrinsic core radio luminosity from the gen- 
eral correlation between the core and the total radio power 
(Giovannini ct al. 1988) we can derive constraints on the 
orientation with respect to the plane of the sky of the 
arcsecond core emission. Objects with a core radio power 



10 



A. p. Schoenmakers et al.: WNB 0313+683: Analysis of a Newly Discovered Giant Radio Galaxy 



Table 2. The measured radio flux densities of the source 
WNB 0313+683 and its components. 



Preq. 


Total 


South Lobe 


North Lobe 


Core 


MHz 


Jy 


Jy 


Jy 


mjy 


38 


15.9 + 1.6 








151 


5.4 + 0.5 ^ 








OA i 


o i9_i_n ^ 


O/I'^-I-OIO^ 
Z.4o zt U.lZ 


n fio -\- n riQ ^ 
u.oy zt u.uo 




343 


3.14 + 0.15 


2.42 + 0.13 


0.72 + 0.03 




1400 


0.91+0.10 " 








1400 


0.82 + 0.03 ^ 


0.63 + 0.02 «5 


0.16 + 0.02 * 


26 + 2 


4850 


0.34 + 0.04 ' 


0.24 + 0.03 j 


0.05 + 0.01 ' 


40 + 3 ^ 


10450 


0.23 + 0.01 ' 


0.12 + 0.01 ' 


0.05 + 0.01 ' 


55 + 3 ' 



Notes: a - 8C survey (Rocs 1990, catalogue revised by Hales 
et al. 1995); b - 6C survey (Hales et al. 1991); c - WENSS ; 
d - Our 92-cm broadband observations; e - NRAO 1400-MHz 
survey flux density (Condon & Broderick 1985) minus NVSS 
background sources flux; f - NVSS; g - NVSS flux density mi- 
nus our 1.4-GHz VLA core flux density; h - Our 1.4-GHz VLA 
observations; i - GB6 survey (Gregory et al. 1996); j - GB6 
minus our 4.85-GHz VLA core flux density; k - Our 4.85-GHz 
VLA observations; 1 - Our 10.4-GHz Effelsberg observations. 
The errors in the 38- MHz and 151-MHz points are estimated to 
be 10%, the errors in all other observations are the quadratic 
sum of an estimated 5% calibration error and a contribution 
due to noise on the map. 



stronger than the expected value are interpreted as galax- 
ies where the core radio emission is Doppler-boosted by a 
relativistic parsec scale jet (see Giovannini et al. 1994 for 
a more detailed discussion) . Comparing the measured core 
emission of WNB 0313+683 with its total radio power we 
can derive an upper limit of 50° for its orientation with 
respect to the plane of the sky. If the radio axis of the 
core emission is the same as that of the large scale radio 
structure, for a projected linear size of this source of 2.0 
Mpc we thus find a real linear size 2.6 Mpc. 



the convention that 5^ oc z/". The maps are shown in Fig. 
9. 

The inverted spectrum core is clearly visible. In both the 
radio lobes there is a gradual steepening when going in- 
wards from the hotspots towards the core. This is best 
visible in the low-frequency spectral index map. Such be- 
haviour is often observed in radio galaxies and can be used 
to estimate the age of the radio source and the advance 
velocity of the lobes (e.g. Leahy et al. 1989; Alexander 
& Leahy 1987). This will be elaborated further in Sect. 
4.2. The central 'bulge' is the region with the steepest 
spectrum, and therefore most probably contains the elec- 
trons with the highest age. The flat-spectrum protrusion 
at the northern edge of the bulge (near Right Ascension 
03'^13™13!52, Declination 68°21'38'.'3) is due to an unre- 
lated radio source visible in the NVSS maps. 
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3.2. Spectral index maps 

Wc have made spectral index maps of WNB 0313+683, 
using the WENSS survey, the NVSS survey, and the Ef- 
felsberg observations. We decided to use the WENSS sur- 
vey map rather than the 92-cm broadband WSRT map 
because it has a much more uniform UV coverage, which 
is reflected in a well-deflned beamshape. This is at the ex- 
pense of having somewhat lower sensitivity as compared 
to the 92-cm broadband observations. 
The NVSS survey maps, which are highest in resolution, 
were convolved to the resolution of cither the WENSS sur- 
vey or the Effelsberg observations. To avoid artefacts, all 
pixels below 3ai in each map were clipped. The spectral 
index a was determined for each individual pixel, using 



Fig. 8. Radio continuum spectrum of WNB 0313+683 for the 
whole source, the southern lobe, the northern lobe and the 
radio core. See Tab. 2 for the flux densities and their references. 



3.3. Polarization properties 

The mere detection of polarized emission in the 92-cm 

broadband WSRT observations already gives an upper 
limit to the Rotation Measure towards WNB 0313+683. 
In this observing mode, each channel has a bandwidth of 
5 MHz, and the highest frequency channel we used is cen- 
tered at 355 MHz. Therefore, if |RM| > 60 rad m^^^ ^j^g 
signal within each channel would have been severely depo- 
larized. Both the 92-cm broadband observations and the 
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Fig. 9. Gr eyscalc plots of the spectral index distribution in WNB 0313+683. a Spectral index between 327 MHz and 1400 MHz, 
with contours from the WENSS map. The greyscale ranges from -0.5 (white) to -1.8 (black), b Spectral index between between 
1400 MHz and 10.45 GHz, with contours from the convolved NVSS map. The greyscale ranges from (white) to -1.1 (black). 



NVSS detect polarized emission over a large fraction of 
the southern radio lobe and the central bulge. This gives 
a total of 6 frequency channels which we can use to find 
and map the RM distribution over the radio source. 



3.3.1. Rotation Measures 

We have used the newstar method (see Appendix A) 
for an analysis of the Rotation Measure distribution over 
the source. Wc used the 92-cm broadband and the NVSS 
data, covering the wavelength range between 20 and 94 cm 
with 6 channels. We have first convolved the NVSS Q and 
U maps to the resolution of the 92-cm WSRT observa- 
tions. In order to uniformly weight each channel, we have 
rescaled the NVSS survey Q and U maps with a factor 
that equalizes the rms noise in the resulting polarized in- 
tensity map to that of the low-frequency WSRT polarized 
intensity maps. 

On the basis of the simulations shown in the Appendix wc 
have decided to remove all pixels in the averaged polarized 
intensity maps with values below IOctp/. Note that even 
with this rejection level we still expect to find points with 



a ~4 rad m~^ offset with respect to the average value. 
A map has been made of the Rotation Measure distribu- 
tion over the source by determining, on a pixel-to-pixel 
base, the RM at which the averaged polarized intensity is 
at its maximum. The RM distribution over the southern 
radio lobe and the central bulge is very smooth. This is 
indicative of a galactic origin. To find the galactic contri- 
bution, we have histogrammed the RMs of the individual 
pixels, and fitted this histogram with a Gaussian. The 
mean of the fit was —10.64 rad m~^ with a standard de- 
viation of 0.23 rad m~^. By subtracting this value from 
the RM-map, we obtain a map of the residual RM, RRM. 
This map is shown in Fig. 10, together with a map of the 
errors in the found RMs. We have defined the error as the 
half-width of the Rotation Measure profile at a polarized 
intensity which lies lapi below the peak of the profile. 
Note that there are indeed a few pixels that have values 
deviating by ~ 4 rad m~^ from their neighbours, as was 
predicted by our simulations. We have not corrected for 
this effect, because the precise value of this deviation is 
not known for each individual pixel, and because only a 
few pixels are involved anyway. 

If the galactic Faraday rotation is uniform over the source. 
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the residual map, as long as its values are significant, in- 
dicates structure which is local to the radio source. We 
find that only a small area towards the central bulge has 
a Rotation Measure which deviates significantly from the 
mean value, with a RRM of —1.9 ± 0.5 rad m~^. 

3.3.2. Depolarization 

We have measured the depolarization distribution be- 
tween the WSRT 92-cm broadband observations and the 
1.4-GHz NVSS radio maps. For the five good channels of 
the 92-cm broadband observations, maps were made of the 
polarized intensity P using P = (Q^ + j/S _ {l.2aQuf'Y''^ ■, 
where aqjj is the average of the noise in the Q and U 
maps. The term containing gqu is a correction for the 
Ricean bias introduced by the quadratic summation of 
maps which each have a Gaussian distribution of noise 
(Wardle & Kronberg 1974). These five maps were then 
averaged into one single map. 

The NVSS /, Q, and U maps were smoothed to the resolu- 
tion of the WSRT 92-cm maps, and maps of P were made. 
All /- and P-maps were set to zero where the flux in the 
92-cm total intensity map did not exceed 3(7/ . Further, 
pixels in the P-maps were blanked when the polarized in- 
tensity did not exceed ^aqu- At each frequency, maps of 
the scalar fractional polarization m! = P/I were made 
(cf. Garrington et al. 1991). From these, the depolariza- 
tion parameter DP}^ = 'm'^43/'m[^QQ was calculated and 
mapped. 

To map the depolarization between the 1.4-GHz NVSS 
and the 10.4-GHz Effelsberg maps, the NVSS maps were 
convolved with a beam of 69" FWHM. The same proce- 
dure as above was followed, but using a cut-off of 3(7/ in 
the NVSS maps. Figure 11 shows the greyscale plots of 
the depolarization between 343 and 1400 MHz, and be- 
tween 1400 and 10450 MHz. At the higher frequencies we 
see only marginal depolarization (DPJ^qq^ 1), but at 
the lower frequencies the depolarization is much stronger. 
Towards the radio core we measure DPgl™ > 1. This 
must be due to the increased contribution of the unpo- 
larized radio core to the total intensity at 1400 MHz. It 
causes a decrease of the fractional polarization and thus 
increases the value of DP^'^^^. 

3.4- Physical properties of the radio source 

We used the measured flux densities at 151 and 327 MHz 
to interpolate the flux density at 178 MHz. Using this flux 
density and the redshift of WNB 0313-1-683, we derive a ra- 
dio power at 178 MHz of Pns = 1.71 x 10^6 W Hz"!. This 
power is comparable to that which distinguishes FRl-type 
from FRE-type radio sources (Fanaroff & Riley 1974), like 
in other Mpc-sized radio sources (e.g. Saripalli et al. 1986). 
We have calculated the equipartition values of the en- 
ergy density and the magnetic field strengths of the entire 
source, the two lobes, and along two slices roughly halfway 



between the radio core and the hotspots in each lobe. 

We use the standard method outlined in Miley (1980), 
but using a conical, rather than cylindrical, geometry for 
each side of the source because of the existence of the 
central bulge. This gives a correction factor of i in the 
volume calculation. We use a base diameter of 840 kpc, 
which is the diameter of the central bulge measured be- 
tween the 6(7/ contours on the 92-cm WSRT map, and a 
length of 2000 kpc, i.e. the projected length. We assume 
an equal amount of energy in relativistic electrons and 
protons (fc = 1), and a filling factor of 1. We use the inte- 
grated fiux densities at 343 MHz of the 92-cm broadband 
observations. The results are shown in Tab. 3. 
From the energy densities the equipartition pressures have 
been calculated assuming that the pressure is dominated 
by relativistic particles. In that case the pressure is just 
^ Ueq- The energy densities and pressures of the lobes are 
best represented by the values resulting from the slices, 
since they do not incorporate the bright and compact 
hotspots. The equipartition energy densities and pressures 
we find are well within the range found in other Mpc-sized 
radio galaxies (Subrahmanyan & Saripalli 1993, Mack et 
al. 1998). We also find low equipartition magnetic field 
strengths, ^ 5 x 10"'' G. Since, at a redshift of 0.0901, 
the equivalent magnetic field strength of IC scattering of 
Microwave Background photons is ^ 3.9/iG, this must be 
the dominant energy loss mechanism for the relativistic 
electrons in the lobes (see Sect. 4.2 for a further discus- 
sion of this topic). 



Table 3. Equipartition parameters of WNB 0313-1-683. All val- 
ues axe calculated assuming equal energies in electrons and 
protons and a filling factor of 1. The energy density is denoted 
by Ueq, the pressure by peq and is just ^Ueq. Beq is the equipar- 
tition magnetic field strength. These three parameters are cal- 
culated for the source as a whole, for the two radio lobes, and 
for two slices perpendicular to the radio axis located halfway 



between the 


core and the hotspot. 






Component 


^eq 


Peg 


B.,q 




10"" erg cm-3 10" 


dyn cm~'^ 


iiG 


Total 


7.7 


2.6 


0.9 


North lobe 


4.5 


1.5 


0.7 


North slice 


1.2 


0.4 


0.4 


South lobe 


11.0 


3.7 


1.1 


South slice 


1.9 


0.6 


0.5 
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Fig. 10. a Greyscale plot of the Residual Rotation Measure. The Residual Rotation Measure is obtained by adding 10.64 rad m~ 

to each pixel. The greyscale ranges from —2.5 to +2.5 rad m~^. b Greyscale plot of the error in the Rotation Measure. The 
greyscale ranges from to 1 rad m^"^. Both maps were obtained using the newstar method with the 92-cm WSRT broadband 
data and the NVSS data combined. 



3.5. Physical properties of the AGN from the optical spec- 
trum 

3.5.1. Extinction determination 

Table 4 lists the measured line strengths of the optical 
spectrum of the host galaxy of WNB 0313+683. Prom the 
ratio of the Ha and H/3 line strengths wc have determined 
the color excess E[B — V) to be 0.90 ± 0.10 mag if it is 
local to the source, or 0.98 ± 0.10 mag if it is caused in 
our own galaxy. We used the interstellar extinction curve 
from Cardelli et al. (1989), and an unreddened Ha/H/3 ra- 
tio of 3.1, as suggested by Halpern & Steiner (1983) and 
Gaskell & Ferland (1984). Since we do not detect a broad 
component of the H/3 line, we only use the narrow line 
components. 

The galactic latitude of WNB 0313+683 is +9?8. From 
the Leiden-Dwingeloo Hl-survey (Hartmann 1994) we 
have extracted the galactic HI column density towards 
the radio source. We find a column density N{HI) of 
2.86 X 10^""^ atoms cm^^, integrated from —450 km s~^ to 
+400 km s-^ Using a conversion factor of N[HI)/E{B — 
V) = 5.6 X 10^^ atoms cm~^ mag~^ (Burstein & Heiles 



1978), we estimate the galactic color excess E{B — V) = 
0.51 mag. We note however that the HI column density 
map of this region is highly structured, so that higher 
density patches may exist within a single beam of the 
Leiden-Dwingeloo survey (~ 30'). Most of the observed 
extinction is thus probably galactic in origin, and we have 
corrected the line fluxes using an E{B — V) of 0.98 mag. 
These corrected line strengths are also printed in Tab. 4. 

3.5.2. Physical properties of the emission-line region 

Figure 12 shows the spectrum in the range surrounding the 
Ha line. Clearly, the line profile of Ha has a broad base. 
The faint broad wings of Ha can be seen to extend from 
roughly 7100A to 7250A, equivalent to ~ 6000 km s'^ 
Ha is the only permitted line with a broad component in 
our spectrum. WNB 0313+683 is not exceptional in this 
case. Some other GRGs have been found to possess broad 
lines, such as WNB 1626+5152 (Rottgering et al. 1996) 
and 0319+411 (de Bruyn 1989). 

Prom the corrected line strengths we can derive the tem- 
perature and density of the narrow- line emitting clouds. 
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Fig. 11. Plots of the depolarization towards WNB 0313+683. a Depolarization between 10.45 GHz and 1400 MHz, with the 
depolarization parameter DP ranging from (white) to 1 (black). Contours from the NVSS maps convolved with the EfFelsberg 
beam (69" FWHM). b Depolarization between 1400 and 343 MHz with DP ranging from (white) to 1 (black). Contours are 
total intensity from the 92-cm broadband observations. 

Prom the ratio of the [OIII] lines, ([OIII]4959 + 5007) / 
[OIII]4363, we find a temperature in the line emitting re- 
gions of 1.61q'3 ^ large error is mainly due to 
the large uncertainty in the line flux of [OIIIJ4363. The 
ratio of the amount of flux in the [SII]6717, 6734 emission 
lines gives an electron density ne of 7^2 x lO'^ cm"'^ in the 
line-emitting regions. 

3.5.3. Optical emission line luminosity and jet power 

Rawlings & Saunders (1991) report a positive correlation 

between the optical emission line luminosity (OELL) of 
radio galaxies and quasars, and their jet power. Although 
the spread in this correlation is large, it holds over several 
orders of magnitude. In their sample they included several 
known GRGs, and they all flt well in this relation. 
In Sect. 4.3 we derive that the total jet power is 4.0 x 
10'''* erg s"-'. The method we use to calculate this is sim- 
ilar to that of Rawlings & Saunders, so the results are 
directly comparable. 

To estimate the total luminosity in optical emission lines 
we also use the same method as Rawlings & Saunders: 
Ltot = 3 X (3 X L[oii]3727 + 1-5 X -^[o//7l4959-i-5007)- From 
Tab. 4 we flnd Ltot = (9.7 ± 0.7) x lO''^ erg s'^. When we 
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Fig. 12. Enlargement of the optical spectrum between 
68OOA and 74OOA, containing the [01], Hq;+[NII] and [SII] 
emission linos. The dashed line roughly indicates the level of 
the continuum around the Ha+[NII] lines 
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Table 4. Measured lines and fluxes in the spectrum of WNB 0313+683. The fluxes are measured using Gaussian flts to the 
observed line profiles. The numbers between the brackets denote the errors in the measured quantities. For the Ha/ [Nil] complex, 
the lines were deblended using three Gaussians; for the [Sll] lines, two Gaussians were used. The powers have been calculated 
using a mean redshift of 0.0901, and assuming isotropic emission. 
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5.2 


(0.4) 


7.1 


[0.4) 


[01] 6364 


6939.1 


0.0904 


9.4 


0.22 


(0.05) 


3.0 


[0.4) 


2.0 


(0.4) 


2.6 


[0.4) 


[Nil] 6548 " 


7139.1 


0.0903 


9.1 


0.82 


(0.2) 


8.7 


[3.0) 


6.3 


(0.8) 


8.7 


[0.4) 


Ha " 


7156.0 


0.0904 


9.1 


3.10 


(0.2) 


32.8 


[3.0) 


23.9 


(0.8) 


32.7 


[0.3) 


[Nil] 6583 " 


7178.3 


0.0904 


8.7 


1.87 


(0.2) 


19.8 


[3.0) 


14.3 


(0.8) 


19.6 


[0.2) 


[Nil] + Hq 








7.9 


(0.4) 


103 


[10) 


60.8 


(3.1) 


84.6 


[4.3) 


[SII]6717 


7323.9 


0.0904 


9.7 


1.23 


(0.04) 


15.9 


[2.0) 


8.8 


(0.3) 


12.3 


[0.4) 


[Sll]6734 


7338.6 


0.0900 


10.8 


1.25 


(0.04) 


16.2 


[2.0) 


8.9 


(0.3) 


12.5 


[0.4) 


[SII] total 








2.51 


(0.17) 


33.4 


[3.0) 


18.2 


(1.2) 


25.3 


[1.7) 



Notes: 

Using an E(B — V) of 0.98 mag. (galactic extinction derived using an Ha/H/3 line-ratio of 3.1). 

Measured using Gaussian fits to the line profiles. 
° The errors do not incorporate the error of 0.10 in the color index; when included, this leads to an additional error of 30% in 
the corrected flux and power. 

Overlaps with the Hgl atmospheric line at 546lA 

Values are calculated using a 3-component Gaussian fit to the (blended) lines. The given Ha flux is for the narrow component 
only. 

^ These values include the broad component of Ha. The total fiux in Ha is 5.17 x lO-i® erg s-i cm-^. 



compare WNB 0313+683 with the sources of Rawlings & 
Saunders, we find that it Hes well away from their correla- 
tion: Its OELL is a factor 10- 20 too high for its jet power 
(or, equivalently, its jet power is a factor 10-20 too low for 
its emission-Une power). We note, however, that when we 
use a galactic color excess E{B — V) of 0.51 mag., as pre- 
dicted by the HI column density from Hartmann (1994), 
WNB 0313+683 agrees much better with the correlation 
of Rawling & Saunders, although the OELL still is on the 
high side. 

4. Analysis of the spectral index distribution 

Jf-.l. Spectral index profiles 

We have integrated the total intensity distribution of 
WNB 0313+683 in boxes perpendicular to the radio axis 
at 327, 1400, and 10450 MHz. The NVSS 1400-MHz data 
were first convolved to the WENSS or EfFelsberg res- 
olution. The width of the slices is half a beamwidth 
(i.e. 30" at 327 MHz, 35" at 10.45 GHz), their length is 



~ 200" (~ 500 kpc). 

Since we are primarily interested in the behaviour of the 

spectral index as a function of distance from the hotspots, 
we have plotted this for each lobe separately. The positions 
of the hotspots have been taken from the high-resolution 
VLA observations. The results arc shown in Fig. 13. 
Between 327 and 1400 MHz, both lobes show a significant 
steepening of their spectrum with increasing distance from 
the hotspot. However, both lobes also show a small but 
significant flattening in the first 100-200 kpc behind the 
hotspots. In the southern lobe, the peak coincides roughly 
with the secondary hotspot that is well visible in the VLA 
maps, but in the northern hotspot there is no such coun- 
terpart. The spectrum near the secondary hotspot in the 
southern lobe is flatter between 1400 and 10450 MHz than 
between 327 and 1400 MHz. Because of the faint high- 
frequency emission of the bridge of the northern lobe a 
spectral index could be accurately measured in this re- 
gion. The strong increase at ~ 800 kpc in the southern 
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lobe and at 1100 kpc in the northern lobe is due to the 
inverted spectrum radio core. 

Jf..2. Expansion velocity, backflows, and the age of 
WNB 0313+683 

The spectral index profiles have been used to estimate the 
age of the radio source WNB 0313+683 and the expansion 
velocities of the two radio lobes. Ageing of a population of 
relativistic electrons results in a steepening of the emitted 
radio spectrum above a certain frequency (e.g. Kardashcv 
1962, Pacholzcyk 1970). The three standard models that 
describe the change in the spectrum are the Kardashev- 
Pacholzcyk (KP), the Jaffc-Pcrola (JP) and the Continu- 
ous Injection (CI) model; see Carilli et al. (1991) for an 
excellent overview of these models. The CI model incorpo- 
rates a continuous injection of 'fresh' electrons, the JP and 
KP models do not and are therefore 'pure' ageing mod- 
els. This, in principle, makes the CI model more suited to 
describe integrated spectra of radio sources, whereas the 
KP and JP models are more suited for analysing spectra of 
source components which do not contain possible sources 
of acceleration. In terms of the spectral shape, the JP and 
the CI model give the most, respectively the least amount 
of steepening. Therefore we decided to use only these two 
models in the following analysis. 

The important parameters that describe the evolution of 
an ensemble of radiating electrons are the time past since 
their last reacceleration, the magnetic field strength, the 
intensity of the microwave background radiation and the 
low-frequency spectral index of the radiation, often called 
the injection spectral index a.inj. Other factors also modify 
the radio spectrum, such as expansion of the radio emit- 
ting plasma. However, our analysis concentrates on the 
bridges of the radio source, where we believe expansion is 
not important any more. 

The equivalent field strength of the microwave background 
radiation (MWBR) is given by 3.24 (1 + z)2 ^G, adopting 
a present-day temperature of 2.726 K for the MWB radi- 
ation. At the redshift of WNB 0313+683 {z = 0.0901) this 
yields Bm = 3.9/iG, which is almost ten times larger than 
the internal equipartition magnetic field strength in the 
radio lobes (see Tab. 3). Therefore IC scattering of the 
MWBR photons must be the dominant energy loss fac- 
tor for the electrons in the lobes of WNB 0313+683. The 
maximum amount of time past since the last acceleration 
of the electrons, which we define as the age of the source, 
is then obtained by assuming that the internal magnetic 
field strength of the source Bg = Bm/V^ (e.g. van der 
Laan & Perola 1969). For WNB 0313+683, an upper Hmit 
on the source age is thus obtained when B^ = 2.1/iG. We 
have used this value in our spectral age calculations. 
To find the advance velocities of the hotspots and the ages 
of the lobes we used the following method. Under the as- 
sumptions given above, we have first calculated a library 
of radio spectra as a function of age of the source and in- 



jection spectral index. Prom these we have derived tables 

of the spectral indices between 325 and 1400 MHz, and 
between 1400 MHz and 10.45 GHz. 

For a fixed advance velocity of the hotspot and amj , each 

position in the radio lobe has a unique age (neglecting 
possible backfiows) and thus a unique spectral index be- 
tween two fixed frequencies. So, using a range of amj and 
advance velocities, we can compare the expected spectral 
index behaviour with our observations. How well the ex- 
pected behaviour fits the observations is measured by cal- 
culating the (reduced) xt- This results in an array of xt 
values. In this array we search for the minimum xti which 
gives the parameters that fit the observations best. The 
errors in the fitted parameters are found by searching the 
la error-ellipse in the x^-aii&y, which for a two-parameter 
problem is outlined by array positions with a value of the 
minimum plus 2.3 (e.g. Wall 1996). This procedure has 
been followed for both JP and CI models. 
Our method differs from the more conventional way 
of measuring advance velocities of hotspots from two- 
frequency radio data (e.g. Myers & Spangler 1985). Usu- 
ally, one assigns a certain amj , often taken as the hotspot 
spectral index, and then calculates the age of each mea- 
sured point in the lobe. When this is plotted against 
hotspot distance a linear fit is made through the measured 
points in order to obtain a velocity. The advantage of our 
method is that both ounj and the hotspot velocity are de- 
rived from the data, and that reliable error estimates for 
both values are obtained. The minimizing of the en- 
sures that the best solution under the given assumptions 
is found. 

One important point to take into consideration is that our 
method only gives physically relevant results if the posi- 
tion of the last acceleration of the electrons is known. This 
is most likely the position near the head of the lobe with 
the flattest spectrum, which usually coincides with the 
primary hotspot. In the southern lobe of WNB 0313+683 
the secondary hotspot has the flattest spectrum between 
327 and 1400 MHz. We have therefore used the position of 
the secondary hotspot as the zero-age point for the radi- 
ating electrons. Also in the northern lobe the region with 
the flattest spectrum does not lie at the head of the lobe. 
Again, we have decided to use the point with the flattest 
spectrum as the zero-age point. 

In the southern lobe, we have fitted both the low- and 
high-frequency spectral index profiles to find out if they 
yield consistent results. In the northern lobe, we only 
used the profile of the spectral index between 327 and 
1400 MHz. The 10.45-GHz flux densities are too low to be 
used for any fitting. However, as a consistency check, we 
have plotted the profiles for the expected spectral index 
distribution between 1400 MHz and 10.45 GHz. 
The best fitting values for amj and the hotspot velocity 
have been presented in Tab. 5. Note that the errors given 
in the table are only the formal errors resulting from the 
fit procedure, and do not take into account possible sys- 
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tematical errors in the data. The spectral index profiles 

resulting from these parameters have been drawn on top 
of the measured spectral index profiles (Fig. 13). 
We find that the best fit of the low-frequency spectral in- 
dex profile in the northern lobe is yielded by a JP model, 
with a lobe advance velocity of 0.027 ± 0.002c. The high- 
frequency profile constructed with this velocity is not in- 
consistent with the measured data points. For the south- 
ern lobe, the JP and CI model fit the spectral index profile 
equally well, both at low and high frequencies. Since the 
low-frequency spectral index measurements have smaller 
errors, we will use those results for an upper and lower 
limit on the hotspot advance velocity in the southern lobe. 
The good fit of the JP model to the northern lobe indicates 
that the radiating particles in that lobe must be relatively 
undisturbed. The southern lobe, with its much higher sur- 
face brightness, brighter and double hotspot, and its cen- 
tral bulge, is much more indicative of a backflow with 
mixing and re-acceleration of the radiating particles. This 
might be why the CI model still provides a good fit to the 
data here. 

If these results have any physical meaning (which, because 

of the large number of assumptions, is questionable; see 
e.g. Eilek 1996), the advance velocities of the hotspots of 
WNB 0313-^683 are between 0.02 and 0.04c. Velocities of 
radio lobes have been determined for a sample of power- 
ful 3C sources by Alexander & Leahy (1987). They find 
velocities in the range of 0.01c - 0.2c, and a weak corre- 
lation of lobe velocity with radio power. WNB 0313-1-683 
agrees quite well with that correlation, since it has a low 
radio power and a low lobe advance velocity. The same 
was found for other GRGs (e.g. Parma et al. 1996, Lacy 
et al. 1993). 

The age of the northern lobe can be estimated by calcu- 
lating the crossing time from the core to the head of the 
lobe (- 1200 kpc) with a velocity of 0.027 ± 0.002c. This 
results in an age of 1.4 ± 0.1 x 10^ yrs. Similarly, for the 
southern lobe, we find an age between 0.7 ± 0.1 x 10^ yrs 
(JP model) and 1.4 ± 0.1 x 10^ yrs (CI model). If back- 
flows in the lobes are important, then this will increase 
the estimated age of the source. 

The assumptions we have made were meant to provide us 
with an upper limit on the source age. We therefore use 
the highest age we find, 1.4 ± 0.1 x 10^ yrs, as an upper 
limit. Ages around 10^ yrs have been found to be fairly 
typical of GRGs (e.g. Mack et al. 1998). 

4-3. The density around the radio lobes (from ram- 
pressure arguments) 

In case that the head of a lobe is ram-pressure confined, 
the density pa of the ambient medium can be found using 
Pa = Ilj/{Afivf^), where Ilj is the thrust of the jet. Ah is 
the area of the bowshock and Vh is the advance velocity of 
the head of the lobe. The thrust Ilj of the jet is given by 
Qjet/vj (e.g. Begelman et al. 1984, their Appendix B.3), 



Table 5. The fitted injection spectral index ainj and the ad- 
vance velocities of the northern and southern radio lobe. For 
the southern lobe, both the low and high frequency spectral 
index profiles were fitted. Velocities are given in units of the 
speed of light c. Also given is the reduced of the fit. The 
last column gives the range of points which were used for the 
fit, in kpc from the hotspot. 



Model ainj v/c xt Fit range 

Southern Lobe 

Low frequency spectral indices 
CI -0.73 ±0.02 0.020 ±0.002 0.65 130 - 700 
JP -0.73 ±0.02 0.040 ±0.003 0.63 

High frequency spectral indices 
01 -0.65 ±0.02 0.024 ±0.005 0.73 100 - 600 
JP -0.66 ±0.02 0.060 ±0.005 0.77 



Northern Lobe 

Low frequency spectral indices 
CI -0.64 ±0.03 0.014 ±0.002 3.78 100 - 1000 
JP -0.62 ±0.03 0.027 ±0.002 0.38 



with Vj the velocity of the material in the jet and Qjet 
the amount of energy delivered by the jet per unit time, 
or the jet power. We can estimate this by dividing the 
total energy content of the radio source by the age of the 
source. 

The energy density of WNB 0313+683 has been measured 
and is given in Tab. 3 as 7.7 x 10^^^ erg cm~'^. To esti- 
mate the volume occupied by the radio lobes, we again 
assume that WNB 0313+683 has a double conical mor- 
phology (because of the central bulge), with a total length 
of 2000 kpc and a base diameter of 840 kpc {6a contours 
in the 92-cm WSRT radiomap) . We thus find a total vol- 
ume of 3.7 X 10^ kpc"^. Therefore the total energy con- 
tent of the radio source, assuming it is in equipartition, is 
8.4 X lO'"^*^ erg. 

An upper limit for the age of the source has been found to 
be 1.4 X 10^ yrs. If the jet power Qj^t has been constant 
during this time, the energy delivered by each of the two 
jets is Qjet ri~^ x 10''"' erg s^^. Here rj is the efficiency 
by which the kinematic energy carried by the jet is con- 
verted into radiation. A conservative choice of rj is given by 
Rawlings & Saunders (1991), who take rj = ^ stating that 
half of the energy is lost to work done by expansion and 
half of the energy is transferred into radiation. Norman 
et al. (1982) have found that rj depends on the density 
contrast between the material in the jet and the exter- 
nal medium. Since the density of the external medium of 
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Fig. 13. Plots showing the spectral index profiles along the radio axis of the two lobes of WNB 0313+683. The two upper 
panels show the spectral index profiles between 327 and 1400 MHz, 092, the lower panels between 1400 MHz and 10.45 GHz, 
a2(f • Because of the lower resolution of the Effelsberg data, the sizes of the bins are somewhat larger in the high-frequency 
plots. Also plotted are the best fitting profiles from the model fits (see text for details). The dashed lines indicate the fits using 
a JP, the dotted lines using a CI model. 
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GRGs is low, the efficiency is probably also low unless the 

density of the material in the jet is low as well. As yet, 
there are no observational constraints on the efficiency, so 
we will use the conservative choice rj = ^ to find a lower 
limit on the jet-power. This results for WNB 0313+683 in 
Qjet ^ 2 X 10''^ erg s~^. The velocity of the material flow- 
ing down the jet we set at c, the speed of light, to obtain 
a lower limit on the thrust. 

The area of the bowshock is difhcult to find since this is 
not directly observable. If we assume that it is not larger 
than the area of the observed hotspots, we can use our 
VLA observations to find an upper limit. To obtain max- 
imal spatial resolution, we made new maps of the VLA 
data using uniform, rather than natural weighting of the 
UV-data. We then fitted both hotspots with a single Gaus- 
sian, using the JMFIT program in the AIPS software pack- 
age. The deconvolved diameter of the southern hotspot 
is 8'.' 5, which translates into a maximum impact area of 
280 kpc^. For the advance velocity of the hotspot we use 
0.03c, which is the mean of the two values from Tab. 5. 
We thus find a lower limit on the density of the external 
medium pa 3.7 x 10~^° g cm~^, or on the particle den- 
sity Ha ^ 1.6 X 10~^ cm~^ (using a mean atomic mass per 
particle of 1.4 amu). 

The hotspot in the northern lobe has a deconvolved width 
of 14", equivalent to an impact area of 780 kpc^. Using 
an advance velocity of 0.027c, we find an ambient den- 
sity pa ^ 1-5 X 10~^° g cm~^ and for the particle density 
Ua 5.8 X 10"'^ cm-3. 

The densities we find are in good agreement with those 
found around other Mpc-sized radio galaxies (e.g. Sari- 
paUi et al. 1986, Parma et al. 1996, Mack et al. 1998). The 
density contrast between the two sides coincides with the 
asymmetry in the hotspot flux: the side with the higher 
luminosity is also the side with the higher external den- 
sity. This agrees with the prediction that the efhciency of 
the acceleration process increases with increasing density 
of the ambient medium (e.g. Norman et al. 1982). 
We have so far assumed that the jet is a straight par- 
ticle stream. In case that it precesses, the area of the 
bowshock in the above analysis should be replaced by 
the much larger cross-section of the entire radio lobe 
(Scheuer 1982). This results in much lower external den- 
sities (~8 X 10~^ cm~^ for the southern lobe). An indica- 
tion for the presence of such a precession is given by the 
observation of two hotspots in the southern lobe. Numer- 
ical simulations of precessing jets reproduce such struc- 
tures very well (e.g. Cox et al. 1991). A valuable, but 
difficult, observation would be to map the jet along the 
entire length of the lobe. 



5. The origin of the observed depolarization 

5.1. The Faraday dispersion 

To find the origin of the observed depolarization, it is con- 
venient to convert the depolarization parameter DP into 
a more physical parameter, the so-called Faraday disper- 
sion, denoted by A. It is defined as the dispersion in the 
Faraday function F(^) = / neB\\dl, the integral over the 
line-of-sight of the product of the electron density and the 
component of the magnetic field along the line-of-sight 
(Burn 1966). By lack of knowledge about the true Fara- 
day function, it is usually assumed that F{(j)) has a Gaus- 
sian distribution, with standard deviation A. It can then 
be shown that, for a source with a redhift of z, the frac- 
tional polarization m at observing wavelength A, m(A) = 
to(0) exp(— 2fc^A^A^(l + z)~^). If A is expressed in units 
of cm-3 pG pc, k = 0.81. In the case of WNB 0313+683 
we thus find A = 1.33(- ln£»P^|§0)V2 cm'^ pG pc. 
Only the southern lobe is well detected in polarized emis- 
sion at 343 and 1400 MHz. To study the behaviour of 
the depolarization in the southern lobe as a function of 
distance from the core, we integrated the total intensity 
and the polarized flux in annuli of 30" width (~ half a 
beamsize), centred on the radio core, at both frequencies. 
We excluded the unrelated radio source at RA 03'^13™50'^, 
Dec. 68°11'50" by blanking it on the maps. Only points 
where the total intensity in the 343-MHz map exceeded 
6 mjy beam^"'^ were used in the integration. For each an- 
nulus, the scalar fractional polarization was calculated at 
each frequency, and from this the depolarization measure 
DPs^f " and the Faraday dispersion A. A plot oiDPj^^§° as 
a function of radius from the core is shown in Fig. 14. For 
reference, the dotted line gives the integrated total inten- 
sity for each anmilus in arbitrary units. Strangely enoiigh, 
the profile of the depolarization parameter DP follows the 
line of total intensity quite well. 

At large radii (J; 700 kpc, near the southern hotspot) 
there is only little depolarization. Towards the core there is 
a small increase, but DP remains more or less constant be- 
tween 200 and 500 kpc. Close to the core, there is again a 
decrease in depolarization, but this is due to the influence 
of the unpolarized core, as was already explained in Sect. 
3.3.2. The plot of the Faraday dispersion A is also shown 
in Fig. 14. Around the hotspot, A 0.85 cm~^ pG pc. To- 
wards the core it increases to a value of ~ 1.5 cm~^ pG pc 
at a radius of 600 kpc, and decreases to a value close to 
1 cm~"^ pG pc at radii between 200 and 500 kpc. Com- 
pared to studies of other radio galaxies (e.g. Johnson et 
al. 1995), the Faraday dispersion towards WNB 0313+683 
is low, even in the innermost 100-200 kpc of the bridges. 
There are three possible regions where the observed de- 
polarization can occur: First, in the halo of our own 
galaxy, second, inside the radio source itself, and third, in 
a large cluster halo surrounding the entire radio source. 
We will discuss all three possibilities for the case of 
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Fig. 14. Radial profiles of (a) the depolarization parameter DPl^'il^, and (b) the therefrom derived Faraday dispersion A 
towards the southern radio lobe. The dotted line in the left plot is the radial profile of the total intensity. The open symbols at 
small radii axe the values after the unpolaxized core has been subtracted from the data. For clarity, they have been shifted by 
a small amount in radius. 



WNB 0313+683. 



5.2. Depolarization by a galactic foreground screen 

Magneto-ionic clouds in the Galactic plane and halo can 
depolarize emission from extragalactic radio sources, and 
cause Faraday effects in the polarization of the syn- 
chrotron radiation from our own Milky way. These ef- 
fects were observed by Wieringa et al. (1993) using the 
WSRT, and show up in their maps as patches of polar- 
ized emission without any related signature in the total 
power maps. They interpreted these observations by as- 
suming that there must be clouds or filaments of magneto- 
ionic plasma between the observer and the smooth galactic 
background emission. Such filaments rotate the polariza- 
tion angle of the incoming radiation, causing observable 
structures in the polarized intensity maps of the WSRT 
interferometer, whereas the total intensity is not changed 
and, due to its large-scale smoothness, remains undetected 
by the WSRT. 

Wieringa et al. (1993) observed that the smallest scale- 
sizes for these filaments were 5'-10', which is smaller than 
our radio source. However, depolarization is caused by gra- 
dients in the Rotation Measure within a single beam. Since 
the filaments are much larger than our beamsize, we can 
exclude the galactic halo as the cause of the observed de- 
polarization. Also, if strong gradients were present in the 
Rotation Measure, they would have shown up in our RM- 
maps. As an extra check we mapped the polarized inten- 
sity over the visible area of these observations. We did not 
detect any large-scale polarized structures, however. 



5.3. Internal depolarization 

Because of the large physical size of WNB 0313+683, there 

is considerable pathlength through the source itself. If 
thermal matter is entrained inside the radio lobes, or if 
the radio lobes themselves are filamented, internal depo- 
larization is likely to occur. A classical test for internal 
depolarization is to observe a breakdown of the linear ro- 
tation angle - relation before a rotation of 90° occurs 
(Laing 1985). In our case, however, most of the rotation is 
of galactic origin, and the values of the residual Rotation 
Measure are not significant (see Sect. 3.3). Prom a Gaus- 
sian fit to the RM distribution over the southern lobe, we 
found a dispersion in the RM distribution of 0.23 rad m~^. 
For this value, the rotation of the polarization angle with 
wavelength goes a.s — 6o = 13.2 [°] , where Oq is the an- 
gle at zero wavelength. For a wavelength of 94 cm, which 
is our lowest frequency channel, the internal rotation of 
the polarization angle is therefore at most ^12°. This is 
not high enough for a measurable deviation from the po- 
larization angle - A^ relation. 

If the depolarization occurs internally, we would also ex- 
pect a correlation between the amount of observed depo- 
larization and the 'depth' of the radio source along the 
line of sight. This requires knowledge of the source struc- 
ture along the line of sight, which we do not have, but we 
can assume cylindrical symmetry around the radio axis. 
Since the radio source is 'fatter' near the center than near 
the hotspots, the amount of depolarization towards the 
center should be larger if it were caused internally. This 
is under the assumption that there arc no large changes 
in the strength or scale-length of the magnetic field or 
the density of thermal matter within the source, which 
can also influence the observed depolarization. Figure 14 
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shows that the amount of depolarization indeed increases 
somewhat at radii below 300 kpc, but this is not highly 
significant. 

Internal depolarization is further indicated by a correla- 
tion between gradients in the Rotation Measure distribu- 
tion and the amount of depolarization (Leahy et al. 1986). 
The RM gradients we find are however small and not sig- 
nificant, so we are not able to test if such a correlation 
exists in our data. We therefore have no conclusive evi- 
dence supporting internal depolarization, but neither can 
we exclude it on the basis of the data available. 

5.4- Depolarization by a, m,agneto-ionic halo 

A third mechanism to d(^polarize radiation is by a large 
halo surrounding the i;uli(j source. A magnetized plasma 
consisting of thermal (10^ — 10* K) electrons and tangled 
magnetic fields (often called a magneto-ionic plasma) de- 
polarizes radiation passing through it, as long as the scale 
length of the tangled B-field is much smaller than the 
observing beam (e.g. Laing 1985, Garrington & Conway 
1991). The beamsize we use (^ 60") translates into a phys- 
ical size of ~ 135 kpc at the distance of WNB 0313-1-683. 
This is much larger than the maximum scale dmax of mag- 
netic field tangling in a cluster halo, if this is set by the 
motion of galaxies through the cluster gas (i.e. 20 kpc; 
Jaffe 1980). Therefore, we most likely do not resolve the 
individual 'Faraday cells' and do not need to correct for 
this effect on the observed Faraday dispersion. 
The Faraday dispersion for radiation passing through a de- 
polarizing plasma is given by A = (n^B|)^/^ {L dmaxY^'^ 
(e.g. Garrington & Conway 1991). Here, L is the path- 
length through the depolarizing plasma in pc. If a large 
magneto-ionic halo is responsible for the observed depo- 
larization, we expect an increase in the observed Faraday 
dispersion at small radii. Typical core radii of the extended 
X-ray emission in both rich clusters and poor groups of 
galaxies are 200-600 kpc (Mulchacy & Zabludoff 1998). 
Radio sources that are at the centers of clusters with 
strong cooling flows and thus high densities in their central 
regions have been found to be strongly depolarized (e.g. 
Taylor & Perley 1993; Taylor et al. 1994) within their core 
radius. However, for WNB 0313+683 there is no strong in- 
crease in A at radius ^ 200 kpc (Fig. 14b). Partly, this is 
caused by the earlier discussed influence of the unpolarized 
radio core. To remove this effect we carefully subtracted 
the radio core from the total intensity data. We then re- 
calculated the fractional polarization in the inner annuli. 
We found that there is still only a small increase in A at 
distances ^ 200 kpc, with A « 1.4 cm~^ /iC pc. There- 
fore, if there is any cluster gas around the host galaxy, it 
is either very centrally condensed like a galaxy halo (e.g. 
Jones & Forman 1984), or it has a very low density, even 
at the center of the cluster. Sensitive X-ray observations 
should be able to disentangle these two possibilities. In 
either case, the heads of the two lobes of WNB 0313+683 



must be currently residing in extremely low-density envi- 
ronments. 

Earlier we estimated the density of the medium surround- 
ing the radio source using the spectral age of the radio 
source and assuming ram-pressure at the head of the lobe 
(Sect. 4.3). Using the formula for A given above, we can 
now try to give an independent estimate. Assuming that 
the radio source resides in a large gaseous halo, the di- 
ameter of this halo must be at least comparable to the 
size of the radio source. We can thus take L 1 Mpc, 
and find {n'^B'^)^/'^ <, 9.9 x 10"^ cm^^ ^Q. To estimate 
the magnetic field strength in the cluster gas, we can use 
the argument that the pressure of the medium outside 
the radio source cannot be larger than that in the radio 
lobes, in order to maintain an equilibrium situation. Us- 
ing equipartition arguments, this is the same as saying 
that the magnetic field in the cluster gas cannot be much 
stronger than the equipartion magnetic field in the radio 
lobe. From Tab. 3, we therefore take a field strength of 
0.5 fiG as an estimate for the mean magnetic field strength 
in the cluster gas. We derive a lower limit for the electron 
density around the southern lobe, averaged along the line 
of sight, of 2.0 X 10~^ cm~^. This is an order of mag- 
nitude higher than the value derived from ram-pressure 
arguments (Sect. 4.3). However, since both methods de- 
pend on many assumptions this is not very surprising. For 
instance, if part of the depolarization occurs internally of 
the source, the density derived above is too high. As we 
mentioned in the last section, we are not able to exclude 
this possibility. 

6. The nucleus: a restarted radio source? 

The properties that distinguish WNB 0313+683 from the 

general population of GRGs are its inverted spectrum ra- 
dio core, and the relatively large power in optical emission- 
lines, as compared to the jet power (see Sect. 3.5.3). It is 
therefore very attractive to believe that these two phe- 
nomena are somehow related. We discuss several scenar- 
ios to explain these observations. However, a few words of 
caution are appropriate: The estimates of the jet power 
are based on many assumptions, such as the validity of 
equipartition and the correctness of our source age esti- 
mate. The power in optical emission-lines, on the other 
hand, depends strongly on the extinction, which we esti- 
mated using the Balmer decrement, but this plays only a 
minor role in the following discussion. 
Relativistic beaming increases an observed flux (e.g. 
Blandford & Konigl 1979) and, depending on which com- 
ponent is beamed, can change the observed spectrum. 
However, for this effect to be important in the core of 
WNB 0313+683, the angle of the radio axis to the line of 
sight must be very small. Also, beamed radio sources, such 
as BL Lacs, are often highly variable. Further, the optical 
spectrum is then expected to show a signiflcant contribu- 
tion of non-thermal (nuclear) continuum, and broad com- 
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ponents of the emitted lines would be much more preva- 
lent. Our present data arc not good enough to exclude 
core flux variability. However, the continuum in the opti- 
cal spectrum is dominated by stellar light, and the broad 
component of the Ha and H/3 lines are certainly not as 
dominant as observed in quasars. This argues against the 
possibility that the radio jet in the core is oriented close 
to the line of sight. 

The inverted spectrum of the core is very suggestive of it 
being a member of the class of Gigahertz Peaked Spec- 
trum (GPS) radio sources, although the peak in the spec- 
trum is at too high a frequency to be classified as such. 
The radio structure of GPS sources has been investigated 
in detail with VLBI (e.g. Snellen 1997), and often show, 
on parsec scales, a morphology very similar to kpc-sized 
radio galaxies. This has led to the suggestion that they 
are very young radio sources, just starting to work their 
way through the ISM of the host galaxy (Mutel & Phillips 
1988). As they grow larger, their peak frequency decreases 
(Snellen 1997). Therefore, the high peak frequency of the 
core of WNB 0313-1-683 is indicative for a very small radio 
source which in the view of Mutel & Phillips must also be 
very young. But how can such a radio source be the core 
of a Mpc-sized and very old radio galaxy? 
This might be explained by a halted and restarted jet 
formation in the nucleus. An estimate of the age of the 
new source can be provided by the observation that the 
hotspots in the southern lobe of WNB 0313-1-683 are still 
very luminous. Once the inflow of jet material stops, the 
hotspots should fade rapidly, on timescales of order 10^ yrs 
(e.g. Reynolds & Begelman 1997). Their presence indi- 
cates that they must still be energized. The light travel 
time from the core to the southern hotspot is almost 
3 Myrs, so the hotspot will still receive jet material for 
at least this amount of time after the central engine has 
stopped producing the jet. The new radio source, which 
we observe as the inverted spectrum core, therefore has 
a time slot of a few Myrs to form itself. Spectral ages of 
other radio sources, such as the well known radio source 
Cygnus A with an estimated age of only 6 Myrs (Carilli et 
al. 1991), indicate that it is well possible to produce even 
large radio sources within a few Myrs. 
Further, the newly formed jet may be much more powerful 
than the jet that formed the old Mpc-sized radio source. 
This could explain the discrepancy between the estimated 
jet power and the emitted power in emission lines: both 
belong to two different epochs in the lifetime of this source, 
and have therefore no relation to each other. 
GPS sources may also be 'frustrated' radio sources, 
embedded in too dense an environment to allow the 
jet to escape (e.g. O'Dea et al. 1991). In the case of 
WNB 0313-1-683, this would indicate that there must have 
been a large change in the gas dynamics of the host galaxy, 
for instance as the result of a merger with a gas-rich 
galaxy. The optical image we have from the DSS (see Fig. 
2b) is of too poor a quality to investigate the morphology 



of the host galaxy for any signs of such a phenomenon. A 
large amount of gas that has been driven into the nuclear 
region provides fuel for the AGN, so this might explain the 
strong optical emission lines as well. What the influence 
on the radio jet would be is however unclear, since there 
is no good theory yet to investigate this. 
To investigate a possible 'rebirth' of the radio core in more 
detail, VLBI observations would be extremely interesting. 
Does the core resemble a small double-lobed radio source? 
Even more interesting would be a search for superluminal 
motion in this source, since this would give an important 
constraint on the orientation of the radio axis. To estimate 
the amount of gas towards the central engine, a study of 
Hl-absorption would be very valuable. 

7. Summary and conclusions 

We have investigated the newly discovered large radio 
source WNB 0313-1-683. This 15' source was discovered in 
the WENSS. The radio core has been identified with a 
weak optical galaxy. It has a redshift of 0.0901 ± 0.0002, 
and thus a projected linear size of 2.0 Mpc. The optical 
spectrum is typical of a narrow-line radio galaxy. The Ha- 
line has a broad (~ 6000 km s~^) component. The color 
index, determined from the (narrow) Ha/H/? line-ratio, is 
large: E{B — V) = 0.98 ± 0.10 mag, where we have as- 
sumed that the extinction is mostly galactic (the galactic 
latitude is -|-9.8°). The HI column density towards the 
host galaxy, measured in the Leiden-Dwingeloo galactic 
HI survey, translates into a color index of 0.51 mag. But 
the large variations observed on scales of 30' in the HI 
survey maps suggest that denser regions may well exist 
within a single beam. 

In addition to the WENSS and the NVSS 1.4-GHz maps, 

we have observed WNB 0313+683 with the 100-m Effels- 
berg telescope at 10.4 GHz, with the WSRT using the 92- 
cm 8-channel broadband system, and with the VLA at 1.4 
and 5 GHz. The VLA and Effelsberg observations show 
that WNB 0313+683 is a core-dominated radio galaxy. 
The radio spectral index of the core is +0.42 ± 0.03, at 
least up to 10.4 GHz, at which frequency it contributes 
~ 25% to the total fiux of the source. The prominence of 
the core, if due to relativistic beaming, suggests that the 
radio axis of the (sub)arcsec core structure is oriented at 
an angle ^ 50° to the plane of the sky. If the radio axis 
of the source as a whole is similar to the radio axis of the 
core, the deprojected linear size would be ^ 2.6 Mpc. 
Spectral index profiles have been made along the radio 
axis. They have been used to estimate the age and ad- 
vance velocities of the two lobes. We find that the source 
has a maximum age of 1.4 ± 0.1 x 10* yrs. The found 
advance velocities of the lobes are (0.027 ± 0.002)c for 
the northern lobe, and ^ (0.03 ± 0.01)c for the southern 
lobe. The ambient densities have been estimated using 
ram-pressure equilibrium at the heads of the jet, and by 
assuming that the impact area is limited to the observed 
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size of the hotspots. We find lower Umits to the parti- 
cle density of 1.6 x 10^*' cm^'' for the material surround- 
ing the southern hotspot and 4.8 x 10"'' cm"*^ around the 
northern hotspot. The contrast in the external density of 
the two lobes is reflected in the lobe luminosities. The 
southern lobe is much brighter (by a factor of ~ 5) and 
appears more confined than the northern one. 
We have studied the Rotation Measure and depolariza- 
tion towards WNB 0313-^683. The Rotation Measures 
have been determined by a new technique, called Rota- 
tion Measure Mapping. We find a very smooth RM dis- 
tribution over the southern lobe, with a mean value of 
— 10.6±0.2 rad m~^. We believe that this is of galactic ori- 
gin. After subtracting this mean value from the observed 
values, we only find significant residual RM values towards 
the central bulge (—1.9 ± 0.5 rad m"^). 
Depolarization towards WNB 0313-1-683 is only important 
at wavelengths above 21cm, since we measure almost no 
depolarization between the Effelsberg and the NVSS data 
but find significant depolarization between the NVSS and 
the 92-cm WSRT data. We have discussed the origin of the 
depolarization and conclude that it could either be caused 
inside the radio source itself, or by a large magneto-ionic 
cluster-like halo surrounding the source. The current data 
do not allow us to distinguish between these two scenar- 
ios. We have investigated the distribution of the depolar- 
ization, and the derived Faraday depth A, towards the 
southern lobe. We find no strong increase of A near the 
core, as would be expected if the host galaxy lies at the 
center of a cluster halo with a dense core, but the val- 
ues are significantly higher than those near the southern 
hotspot. Overall, A is close to 1.0 cm"'^ /iG pc. 
In the case that a large halo surrounding the radio source 
is responsible for the depolarization, we estimate a lower 
limit for the electron density, averaged along the line of 
sight towards the southern lobe, of 2 x 10"^ cm""^. This 
is an order of magnitude above the values we find in the 
spectral index analysis. Part of this might be caused by 
a contribution from internal depolarization. Our data do 
not exclude this possibility. 

In many respects, WNB 03134-683 is like any other GRG. 
Its size is close to the median value of the 'classical' 
GRGs (~ 2 Mpc), the density of the ambient medium is 
low, of the order of 10~^ cm"-^, the spectral age is large 
(~ 10^ yrs), there is only little depolarization, and the Ro- 
tation Measures, corrected for the galactic contribution, 
are low as well 2 rad m~^). 

The inverted spectrum radio core and the very strong opti- 
cal emission lines as compared to the estimated jet power, 
suggest an interrupted radio activity of the AGN. VLBI 
observations are necessary to investigate the structure of 
the radio core in more detail. 

Acknowledgements. KHM was supported by the Deutsche 
Forschungsgemeinschaft, grant KL533/4-2 and by European 
Commission, TMR Programme, Research Network Contract 



ERBFMRXCT96-0034 "CERES". LL acknowledges Bill Cot- 
ton for useful suggestions on the use of NVSS data to study gi- 
ant radio sources. We would like to thank the referee, U. Klein, 
for many helpful suggestions which considerably improved the 
paper. NRAO is operated by Associated Universities, Inc., un- 
der cooperative agreement with the National Science Founda- 
tion. The INT is operated on the island of La Palma by the 
Isaac Newton Group in the Spanish Observatorio del Roque 
de los Muchachos of the Instituto de Astrofisica dc Canarias. 
The Westerbork Synthesis Radio Telescope (WSRT) is oper- 
ated by the Netherlands Foundation for Research in Astronomy 
(NFRA) with financial support of the Netherlands Organiza- 
tion for Scientific Research (NWO). This research has made 
use of the NASA/IPAC Extragalactic Database (NED) which 
is operated by the Jet Propulsion Laboratory, California Insti- 
tute of Technology, under contract with the National Aeronau- 
tics and Space Administration. 



A. A new method for finding Rotation Measures 

The 'classical' method for finding Rotation Measures is to 
determine the position angle of the polarized emission in 
each individual channel, and to plot these as a function of 
the square of the wavelength. The RM is then the slope 
of the line that fits these datapoints best. This method 
has proved itself extensively, but it provides little insight 
in the spatial distribution of the Rotation Measure when 
studying extended structures. 

For this purpose, a new routine has been developed in the 
NFRA NEWSTAR reduction and analysis package, called 
Rotation Measure Mapping. This routine combines the 
maps of the Stokes Q and U parameter of each separate 
frequency channel in the following way. Using an 'a priori' 
guess of the RM, the Q and U vectors in each channel are 
rotated 'backwards' over the angle that they would have 
been rotated over by the assumed RM, relative to the fre- 
quency of the highest frequency channel. Then, the Q (and 
U) maps of all channels are vector summed. Only if the 
chosen RM is close to the real RM, the rotated Q (and U) 
vectors in each channel will be the same, and the length of 
the summed vectors will be maximal. Using these vector 
summed Q and U vectors, which we will refer to as Qs 
and Us, respectively, we can calculate the summed polar- 
ized intensity Pig = ^/Qf+Uj, which will also be at its 
maximum when the correct RM is chosen. Using a range 
of possible RMs, the RM for which the summed polarized 
intensity is maximal, must be closest to the real value. 
The main advantage of the newstar method is that maps 
of summed polarized intensity can be made for a wide 
range of RMs. These can be displayed in a movie form, so 
that for extended sources, like WNB 0313-1-683, one can 
look for systematic changes in the RM over the source. 
But it can also be used for a more strict analysis of the 
data. 

In the classical method, the spectral index and depolar- 
ization of the emission, when using widely separated fre- 
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Fig. 15. Plots showing the result of the simulations of the NEWSTAR method and the 'classical' method for finding Rotation 
Measures (see text for details). Each plot is a realisation of 500 points with random signal-to-noise values. We used a RM of 
10 rad m~^ as input. Both methods give very similar results. 



quencies, will affect the error in the measured polarization 
angle, and thus the accuracy of the result. In what manner 
exactly strongly depends on the nature of the source ob- 
served. The intrinsic degree of polarization for synchrotron 
radiation with a powcrlaw spectrum is not frequency de- 
pendent. Therefore, in radio sources with negative spec- 
tral slopes, e.g. Si, oc u^^-^, the drop in total intensity at 
a higher frequency can, at least partly, be counteracted by 
the lower amount of depolarization at that frequency. 
In order to test the accuracy and reliability of the new- 
star method, we have performed computer simulations, 
in which we analysed the same artificial input data set 
with the NEWSTAR and the classical method. As input we 
provided 500 random Q and U intensities, which we as- 
sumed to have a frequency of 1400 MHz (cf. the NVSS 
survey). For each datapoint, we then determined the po- 
larization angle at 1400 MHz and calculated the polariza- 
tion angles at the other frequencies, using a fixed Rotation 
Measure of 10 rad m~^. Prom this, the Q and U intensities 
at each frequency have been determined. Then noise was 
added using a Gaussian noise distribution with a standard 
deviation equal to the rms noise of our real Q and U data. 
We find that both methods recover the input RM very 
well, as long as the polarized intensity ^ IQapi, where 
a PI is the noise in the averaged polarized intensity map 
corrected for Ricean bias (see Fig. 15). For polarized in- 
tensities ^ 8<Jpi, the input RM is not reliably recovered. 
Especially note the two 'sidelobes' in the RM distribution, 
which arc separated by ~ 4 rad from the input RM. 
They are the result of the rnr rad ambiguity of the deter- 
mination of the polarization angle in the high frequency 
(NVSS) data. Due to the added noise, chosing a value 
which is 180° higher or lower sometimes gives a better fit 
through the datapoints. We find that 10% of the data- 



points in our simulations which have polarized intensities 
between 8 and 30 api lie in these two sidelobes. 
We have used this method to correct for the observed po- 
sition angles of the E-^eld in the 1400-MHz NVSS ob- 
servations. The result is shown in Fig. 16. This should 
be compared with the orientation of the magnetic field 
that is observed at 10.45 GHz with the Effelsberg tele- 
scope (i.e. the vectors plotted in Fig. 4a rotated by 90° ). 
The ~4 rad m~^ ambiguity causes the small wiggle in the 
B-field vectors in the southern lobe of WNB 0313-1-683. 
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Fig. 16. Contour plot of the WNB 0313+683 at 1400 MHz 
from the NVSS survey (see also Fig. 1). The orientation of 
the vectors represent the orientation of the 'true' B-field after 
rotating the observed _B-ficld with the amount of rotation ex- 
pected from the RM-analysis described in this Appendix. The 
position angles are almost identical to those inferred from the 
high-frequency Effelsberg observations. The apparent 'wiggle' 
is caused by the ~ 4 rad m~^ uncertainty that is described in 
the text. 
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